Evidence for a Spin Phase Transition at v = in Bilayer Graphene 
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The most celebrated property of the quantum spin Hall effect is the presence of spin-polarized 
counter-propagating edge states^™ This novel edge state configuration has also been predicted to 
occur in graphene when spin-split electron- and hole-like Landau levels are forced to cross at the 
edge of the sample 4 ^. In particular, a quantum spin Hall analogue has been predicted at v — in 
bilayer graphene if the ground state is a spin ferromagnelP^. Previous studies have demonstrated 
that the bilayer v — state is an insulator in a perpendicular magnetic fiel d 10 " 16 *, though the exact 
nature of this state has not been identified. Here we present measurements of the v — state in a 
dual-gated bilayer graphene device in tilted magnetic field. The application of an in-plane magnetic 
field and perpendicular electric field allows us to map out a full phase diagram of the v — state 
as a function of experimentally tunable parameters. At large in-plane magnetic field we observe 
a quantum phase transition to a metallic state with conductance of order 4e 2 //i, consistent with 
predictions for the ferromagnet. 



Under a strong perpendicular magnetic field, bilayer 
graphene (BLG) develops a v = quantum Hall (QH) 
state at the charge neutrality point (CNP) which displays 
anomalous insulating behavior^"^. Transport studies 
in a dual-gated geometr^HHI indicate that this gapped 
state results from an interaction-driven spontaneous sym- 
metry breaking in the valley-spin space 17 . However, the 
exact order of the resulting ground state remains contro- 
versial. In high-mobility suspended BLG devices, a bro- 
ken symmetry state at charge neutrality has also been ob- 
served at zero magnet ic fie ld whose nature ha s been un- 
der intense theoreticaPE^U and experimenta l 12 * 14 ^ 16 * 23 * 
investigation. Some experiments indicate ther e may be a 
connection between these two insulating states^^. 

The v = state in BLG occurs at half filling of the 
lowest Landau level. Although the BLG zero-energy 
Landau level has an additional orbital degeneracy stem- 
ming from two accidentally degenerate magnetic oscil- 
lator wavef unctions 2 -^, minimization of exchange energy 
favors singlet pairs in the orbital space at all even fill- 
ing factors^. This leaves an approximate SU(4) spin 
and pseudospin symmetry for the remaining ordering of 
the ground state. Comparable ground state competi- 
tion has also been studied in double quantum wells in 
GaAs 2-dimensional electron systems (2DESs) 26 . Quan- 
tum phase transitions, in particular from canted antifer- 
romagnetic to ferromagnetic ordering, can occur in those 
systems as an in-pl ane m agnetic field is applied (Zeeman 
energy is increased ) 27 * 28 *. 

Our experiment is carried out in BLG samples with top 
and bottom gates in which thin single crystal hBN serves 
as a high quality dielectric on both sides (Fig. 1 bottom 
inset). By controlling top gate voltage (Vtg) and bot- 
tom gate voltage (Vbg) we can adjust carrier density n 
and perpendicular electric displacement field D indepen- 
dently. Additionally, tilting the sample in the magnetic 
field allows us to independently control the Coulomb en- 
ergy in a Landau level, e 2 where the magnetic length 



£b is determined by the perpendicular magnetic field B±, 
and Zeeman energy (determined by the total magnetic 
field B tot ). Since in BLG B tot and D separately couple 
to spin and pseudospin, respectively, our experimental 
setup allows us to tune anisotropics and characterize the 
broken- symmetry QH states in the approximate SU(4) 
spin-pseudospin space (see supplementary information). 

Fig. 1 shows the resistance measured as a function of 
the two gate voltages. At zero magnetic field, the resis- 
tance at the CNP depends monotonically o n \D \ due to 
the opening of a single particle energy gap! 29 * 30 *. With 
increasing applied magnetic field (Fig. 1 top inset) our 
device shows the same behavior reported elsewhere^^yj, 
with a layer-polarized (LP) insulator (which we denote 
phase (i)) observed at large applied displacement field, 
separated from a low displacement field and high mag- 
netic field quantum Hall insulator (phase (ii)) by a well 
defined boundary of finite conductance. For the remain- 
der of the paper we focus on the behavior of phase (ii) 
as we vary the relative Zeeman versus Coulomb energy 
scales in a rotated field geometry. 

Fig. 2a shows a series of conductance maps in the 
vicinity of the CNP. Each pane corresponds to a fixed 
perpendicular applied magnetic field of 1.75 T, but with 
varying total field. At a total field of 1.75 T (that is, 
B± — B tot ), the distinct insulating regimes at low and 
high \D\ are clearly visible, as discussed in Fig. 1. At 
low \D\ the four-terminal conductance increases dramat- 
ically with increasing total magnetic field B to t, reaching 
a value of around 3 — 4 e 2 //i, after which it becomes 
relatively insensitive to either displacement field or mag- 
netic field (Fig. 2b). This rapidly increasing behavior 
of conductance at fixed B± is consistent with the grad- 
ual closing of a gap in phase (ii) as B tot increases. At 
fixed Coulomb energy, this decreasing gap with increas- 
ing Zeeman energy suggests that phase (ii) corresponds 
to a spin-unpolarized ground stat d 8 * 10 * ^ ^. In addition, 
the saturation of conductance at a metallic value around 
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FIG. 1. Device summary. Plot of R xx versus top gate 
voltage Vtg and bottom gate voltage Vbg at B± = B to t = 
5 T. Data taken at 1.5 K. Dashed line indicates the CNP. 
Top right inset shows four-terminal conductance at the CNP 
plotted against B± and D. Dashed line in inset corresponds 
to dashed line in main plot. Bottom left inset shows an optical 
image of the device used in this experiment with false coloring 
of the graphene (pink). 



~ 4 e 2 /h at B tot > 15 T strongly suggests a transition 
into a new phase with a fully closed transport gap, which 
we label phase (iii). 

The transition between insulating and metallic behav- 
iors can further be confirmed by investigating the tem- 
perature dependence of the conductance. Fig. 2 c and d 
plot four-terminal conductance G at the CNP against D 
at fixed B± = 2.5 T with B tot = 2.5 T and 13.1 T, respec- 
tively. Temperature T varies between 3-18 K. At small 
\D\, G exhibits weakly insulating behavior for B to t — B±, 
while at B tot ^> B±, metallic temperature dependence 
(increasing G with decreasing T) is observed. Fig. 2e fur- 
ther contrasts the metallic versus insulating classification 
in different phases by plotting conductance at the CNP 
against temperature for both large D and D = 0. G(T) 
shows insulating behavior for phase (ii) (B± « B tot ) and 
metallic behavior for phase (iii) (B± <C B tot ). At in- 
termediate values of B to t the conductance is relatively 
temperature insensitive, suggesting a gradual crossover 
between the two phases. Interestingly, the cusps of high 
conductance that mark the transition between the LP 
and non-LP phase show metallic temperature depen- 
dence even in a fully perpendicular field. At high D 
(i.e. in phase (i)), the conductance remains relatively 
insensitive to the applied in-plane magnetic field. This is 
expected for a LP v — QH state, where the energy gap 
is determined by the applied displacement field. 

A general theoretical analysis 8 proposes four possible 
ground states for v = in BLG: partially layer-polarized 
(PLP), fully layer-polarized (FLP), canted antiferromag- 



FIG. 2. Tilted field transport, a) 1/R XX plotted against 
carrier density, displacement field, and total magnetic field 
with B± fixed at 1.75 T. Data taken at 350 mK. b) Four- 
terminal conductance at the CNP (G — 1/R XX ) plotted 
against displacement field D for a variety of total magnetic 
fields, B± fixed at 1.75 T. Dashed line indicates the zero dis- 
placement field points plotted in c) . Solid line marks the the- 
oretically expected 4e 2 /h conductance for the F phase. Data 
taken at 350 mK. c-d) Temperature dependent four-terminal 
conductance at the CNP plotted against displacement field, 
with B± = 2.5 T, B to t =2.5 T (c), 13.1 T (d). Colored curves 
are taken at 3 K, black curves at 18 K, and intermediate 
traces are spaced by 1.5 K. Dashed (solid) vertical lines at 
a fixed displacement field correspond to dashed (solid) lines 
in e). e) Four-terminal conductance at the CNP for differ- 
ent total magnetic fields plotted against temperature, all at 
B± — 2.5 T. Solid lines demonstrate how behavior changes as 
the total field is increased for zero displacement field. Dashed 
lines demonstrate similar behavior for all total magnetic fields 
at D = 0.085 V/nm. 



net (CAF), or ferromagnet (F). Considering that layer 
pseudospin is synonymous with sublattice and valley in 
the lowest Landau level 24 , these four different scenarios 
can be visualized in terms of spin and sublattice of BLG 
(Fig. 3a). Which state the system favors depends on 
both anisotropics in the electron-electron and electron- 
phonon interactions, and external symmetry-breaking 
terms like Zeeman splitting ez and interlayer potential 
energy ey. The PLP/FLP phases are favored in the large 
€y limit, and the F phase is favored in the large ez limit. 
Unlike in GaAs double wells, these ground states may 
be accompanied by dramatically different transport sig- 
natures. The BLG ground state configurations discussed 
above are all predicted to be insulators except for the 
F phase. For the F phase, the highest filled band ac- 
quires an electron-like dispersion at the edge while the 
lowest unoccupied band acquires a hole-like dispersion 
(Fig. 3b)R This forces a level crossing at the edge which 
should result in metallic conductance of 4e 2 /h du e to 
spin-polarized counter-propagating edge states^EHl. i n _ 
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FIG. 3. Candidate states, a) Spin and pseudospin configu- 
rations in the BLG v — state. Cartoons represent the layer 
polarized (LP) (i), canted antiferromagnet (CAF) (ii), and 
ferromagnet (F) (iii) phases, b) Diagram of predicted bulk 
and edge energy levels for the CAF (ii) and F (iii) phases. 



terpreting our results in this framework, we can identify 
phase (i) with either the PLP or FLP state, phase (ii) 
with the CAF state, and phase (iii) with the F state. 

We now focus on the phase diagram of the different 
ground states classified above. The phase transitions be- 
tween the LP, CAF, and F phases can be mapped out by 
locating the boundary between the metallic (G > e 2 /h) 
and insulating (G <C e 2 /h) regimes as a function of D 
and B tot with fixed B±. For a given B to t-> we can de- 
fine the critical displacement field for transition D c as 
the D value that gives maximum change of dG/dD in 
G(D,Btot) at fixed B±. Fig. 4a shows the values of D c 
and B to t extracted from Fig. 2b (B± = 1.75 T). In this 
D-Btot ma P three different phases can be identified, cor- 
responding to the LP, CAF, and F phases. Note that 
the transition between the CAF and LP phases is sharp, 
as the conductance changes by an order of magnitude 
over a D variation AD < 0.02 V/nm. By contrast the 
transition from CAF to F is gradual. This observation is 
consistent with the theoretical expectation that the LP 
to CAF or F transition is first order while the CAF to F 
transition is second ordepH 

By identifying characteristic energy scales for the ex- 
ternal symmetry breaking terms, we can comment quan- 
titatively on the phase boundaries for a fixed B±. We 
identify Zeeman energy ez = g^BB to t, where g and 
Hb are the electron spin g-factor and the Bohr mag- 
neton respectively, and the interlayer potential energy 
ey = Dde/2, where d is the interlayer distance in BLG. 
For B± = 1.75 T, at low B tot the CAF-LP transition 
occurs at ey « 5 meV, and at ey = the CAF-F transi- 
tion occurs at ez ~ 0.8 meV. The F-LP transition occurs 
along a line whose slope ey /ez ~ 2. The fact that ey /ez 
at the F-LP phase boundary is of order unity is consistent 
with the increasing stability of a spin-polarized state with 
increasing magnetic field, providing further confirmation 
that the conducting v = phase is a ferromagnet. This 
observation is in qualitative agreement with theoretical 



predictions^. 

Since at the CNP and D = 0, Bj_ is the only phys- 
ical quantity to control the strength of e-e interactions, 
the phase boundary between the CAF and F phases can 
be further investigated by observing G as a function of 
B to t with different values of fixed B± (see supplementary 
information). Fig. 4b shows G(B to t) for three different 
fixed values of B± at D = 0. Generally, the behavior 
of G can be characterized by two different regimes. For 
B to t ~ B± (i-e. the CAF phase), G increases exponen- 
tially as Btot increases. For B to t ^> B± (i.e. the F phase), 
G is insensitive to B to t, saturating close to 4e 2 /h. We de- 
fine the critical magnetic field Bt ot which separates these 
two regimes by the intersection of guide lines for the ex- 
ponential and saturated regimes, as marked by the arrows 
in Fig. 3b. The identification of B^ ot at different B± al- 
lows us to map the phase boundary between CAF and F 
as a function of B±. The inset of Fig. 4b displays B^ ot 
as a function of B±. Although there are only a few data 
points available to indicate the phase boundary, the gen- 
eral trend of the CAF-F phase transition can be inferred 
from them. 

Finally, we remark on the similarity between the F 
phase of BLG at v = and the quantum spin Hall ef- 
fect (QSHE) in HgTcPSl. In the QSH phase quantized 
conduction is due to spin-polarized counter-propagating 
edge states, while the bulk remains incompressible. An 
analogous sce nario is expected for BLG in the F quan- 
tum Hall stat e 4 * 5 * 7 !, except that BLG carriers an addi- 
tional orbital degeneracy. The F state in BL graphene 
is therefore expected to exhibit a four terminal conduc- 
tance of 4e 2 //i, compared with 2e 2 /h for HgTe. Simi- 
lar to reported measurements of Hg Te 3 * 31 !, our devices 
show imperfect quantization with measured conductance 
values less than Ae 2 /h in the F phase. The reduced 
conductance may be due to backscattering in the edge 
state or through remaining impurity states in the bulk. 
Unambiguous experimental support for the existence of 
spin-polarized counter-propagating edge states should be 
provided by well- quantized nonlocal measurements 3 -^ (see 
supplementary information). 

In conclusion, we have measured dual-gated graphene 
bilayers at the v — state. In a perpendicular mag- 
netic field, we observe a transition between two incom- 
pressible states at a finite displacement field, indicating 
that the D = state is not layer polarized. At low dis- 
placement fields, we observe four-terminal conductance 
increase and then saturate as in-plane magnetic field in- 
creases, showing that the D = 0, B to t — B± state is 
also not spin polarized. This is consistent with a quan- 
tum phase transition driven by in-plane magnetic field 
between the canted antiferromagnet quantum Hall state 
and the ferromagnet. The near quantization of the 4- 
terminal conductance to 4e 2 /h agrees with predictions 
that the ferromagnetic quantum Hall ground state car- 
ries spin-polarized counter-propagating edge states. 
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FIG. 4. Phase Transitions, a) Different phases of the v — state mapped as a function of D and B to t at B± = 1.75 T. At low 
Btot, the phase boundary between LP and CAF is determined by the conductance spike separating the two insulating phases. 
At large Btot, the phase boundary between LP and F is determined by the point where conductance begins to exponentially 
decrease with displacement field. The gradual transition between CAF and F with increasing Btot is marked by the gradual 
color change, with F stabilizing at high Btot. b) Four-terminal conductance at the CNP with D = 0, plotted against B to t for 
three different values of B±. Dashed lines indicate regions of exponential increase of conductance and the dotted line marks 
saturated conductance. The critical magnetic field B% ot for the phase transition is marked by arrows. Inset shows Bt ot as a 
function of B± separating the CAF and F phases in a B to t versus B± parameter space. 
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excluded. Four-terminal conductance at the CNP is ob- 
tained by taking G = 1/R XX . 
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